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1.  Research  on  Numerical  Methods  for  Singularly  Perturbed  Differential 
Equations 

During  the  period  of  this  contract,  we  developed  and  applied  numerical  methods 
for  singularly  perturbed  two-point  boundary  value  problems  and  time-dependent  par¬ 
tial  differential  equations.  Our  progress  and  publications  in  these  areas  are  described 
below . 

1.1.  Two  Point  Boundary  Value  Problems 

Our  research  during  this  period  has  largely  centered  on  analyzing  the  asymp¬ 
totic  behavior  of  solutions  to  scalar  problems  of  the  form 

(  y  "  +  f  {y  )y'  +  g  (y)  ^  o.  o  ^  x  ^  r 

with  1/(0)  =.4  and  y(r)=B  prescribed.  Such  problems  and  their  generalizations 
model  a  variety  of  nonlinear  oscillators  and  other  physical  systems  involving  narrow' 
boundary  and  interior  layers  where  solutions  change  rapidly.  A  most  fascinating 
e.xample  involves  /  {y)  =  y(l-y~)  and  =  -  j/ •  The  unique  solution  differs  in 

twenty-two  different  region^;  of  the  .A-B  plane,  having  up  to  four  different  limiting 
solutions  in  each  (depending  on  the  interval  length  c).  .\way  from  regions  of  nonuni- 
forrn  convergence,  the  limiting  solution  I  qIj')  will  satisfy  one  of  the  five  limits; 

(1-  )\-  )  r„=  :  <1-  r,f  )  V„=l.ro(0)=-l;  ll- i'oMTo=l.T()(<' )=B: 

(1-  l'„-  I  }■„=  1.  l  (,(c  )  =  1:  or  TolJ"  1  =  0. 

This  behavior  and  the  location  of  jump  jioinf'  ran  h('  exi)lained  by  analyzing  trajec¬ 
tories  for  the  corresponding  fast-slow  system  in  the  y-z  phase  plane,  where  z  is  the 
Lienard  variable  ;  =  <  y  +F  (y  ].  F  (y  )  =  J  f  (ti  )<L<  .  Mon'  romj)licated  behavior  will 
be  obtained  for  systems  with  more  rest  i)oints  and  when  the  "characteristic  curve" 
:  =  F  {y  )  has  less  symmetry  and  more  extreme  jwints.  Our  rc'sults  may  be  typical  of 
nonlinear  problems  which  arise  as  react i«)n-diffusion  systems  or  as  otlu'r  problems 
with  slow/fast  dynamics.  They  would  be  nearly  imj)ossil)le  to  comprehend  and  dis¬ 
cover  computationally,  without  some  prior  knowledge  of  jxissibh'  analytical  behavior. 
This  work  is  being  carried  out  with  gra<luai<'  student  .John  .\llen.  It  has  benefitt'd 
from  graphical  ODE  routines  j)roduce<l  by  (i.  Odell.  Further  comjtutational  efforts 
will  benefit  from  use  of  adafttive  immerical  methods,  as  have  been  developed  by 
Flaherty.  Brown  and  Lorenz,  and  others.  :ind.  ultimatt'ly.  l)y  (h'vising  domain  decom¬ 
position  techni<iues  relate<l  to  identification  of  asynijitot ic  scales. 

Vasil'eva  and  Butuzov,  and  Flaherty  :ind  O'Malley  have  independently 
developed  analyses  for  singular-singular  perttirt>;ition  {troblems.  Thest-  problems 
occur  when  the  limiting  <’<iuation  has  a  famih  of  solutions.  Tlu'y  arist'  in 
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applications  ranging  from  chemical  kinetics  to  circuit  theory.  Recently,  proposed 
simulation  methods  for  rapidly  oscillating  controls  made  clear  the  need  for  a  more 
comprehensive  theory.  This  is  the  focus  of  thesis  work  by  Gu  Zhong-Mei  and  of  a 
collaboration  with  Dr.  Nikolai  Nefedov,  a  visiting  faculty  member  from  Moscow 
State  University.  At  this  time,  it  appears  that  the  Flaherty-O’Malley  procedure 
(which  provides  a  differential  equation  for  the  limiting  solution)  is  most  practical. 
However,  more  practical  means  of  decoupling  fast  and  slow’  parts  of  solutions  and  of 
examining  appropriate  stability  manifolds  for  stretched  systems  are  both  very  much 
needed.  Results  obtained  will  form  a  basis  for  regularization  methods  to  solve 
differential-  algebraic  systems,  as  well  as  the  original  singular  problems. 

Other  efforts  by  Professor  O'Malley  include  continuing  interactions  with  G. 
•Soderlind  of  Stockholm's  Royal  Institute  of  Technology.  R.  Mattheij  of  the  Catholic 
University  of  Nijmegen,  and  w’ith  colleagues  J.  Cole.  J.  Flaherty,  H.  Kurland,  and  M. 
Slemrod.  among  others.  \  fairly  extensive  set  of  class  notes  on  singular  perturbation 
methods  was  written  last  semester.  It  may  soon  form  the  basis  of  a  research  level 
monograph  or  graduate  textbook. 

1.2  Initial  Boundary  Value  Problems  for  Partial  Differential  Equations 

Professor  Flaherty  and  several  colleagues  and  graduate  students  have  been  con¬ 
ducting  research  on  adaptive  finite  difference  and  finite  element  methods  for  partial 
differential  e(]uations, 

Drew  and  Flaherty  j2''  created  a  model  for  shear  hand  instabilities  that  involves 
the  raj)i(l  shearing  of  a  slab  of  a  visco-elastic  material  with  small  viscosity.  They 
used  a  moving  mesh  finit<'  element  code  to  solve  their  model  and  illustrated  several 
])ossii)le  iiK'chanisnis  for  shear  hand  formations.  In  a  somewhat  related  study.  .Slein- 
rod  and  Flaherty  [l2]  considered  appropriate  finite  difference  schemes  for  phase  tran¬ 
sit  ion>  involving  van  der  Waals  fluids. 

Moving  mesh  methods  were  studied  and  the  stability  of  several  schemes  that  are 
based  on  equidistributing  a  local  error  indicator  was  characterized.  Our  findings, 
reported  in  jlll.  show  that  many  popular  mesh  moving  strategies  can  be  unstable  in 
certain  situations. 

Moving  grid  methods  are  effective  at  reducing  disjiersive  errors  in  the  vicinity  of 
wave  fronts:  however,  no  method  that  uses  a  fixed  numln-r  of  computational  cells  can 
be  tised  to  compute  a  solution  to  a  prescribf'd  level  of  accuracy.  For  this  reason,  we 
have  been  studying  and  developing  techni(}ues  that  perform  exjilicit  error  estimations 

*  Sen  lisl  of  piil)li(-al  ions  and  maniisrnpts  in  Section 
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that  are  used  to  add  and  delete  elements  as  the  temporal  integration  progresses.  In 
this  spirit,  Adjerid  and  Flaherty  [13)  developed  a  finite  element  method  of  lines 
where  the  mesh  "lines"  are  adaptively  moved  to  minimize  artificial  diffusion,  but 
where  elements  are  added  and  deleted  during  the  integration.  Piecewise  linear  finite 
elements  are  used  to  calculate  the  numerical  solution  and  quadratic  elements  are 
used  to  estimate  its  error.  Mesh  motion  is  accomplished  using  one  of  the  stable 
schemes  studied  in  [llj.  Ordinary  differential  equations  are  obtained  for  the  finite 
element  solution,  the  error  estimate,  and  the  mesh  coordinates  and  these  may  be 
solved  using  existing  software  for  stiff  systems.  VVe  have  applied  this  procedure  to 
several  problems,  including  some  difficult  combustion  studies,  and  found  some  very 
encouraging  results. 

The  elementary  global  refinement  strategy  that  was  used  in  jl3]  could  be  ineffi¬ 
cient  in  many  situations.  Additionally,  the  mesh  moving  algorithm  of  [13]  required 
users  to  specify  a  parameter  that  controlled  relaxation  times  to  an  "ideal"  mesh  that 
equidistributes  local  error  estimates.  We  have  remedied  both  of  these  deficiencies 
and  are  now  using  an  efficient  local  refinement  strategy  and  an  adaptive  procedure 
for  selecting  the  mesh  relaxation  parameter.  The  details  of  these  algorithms  and  the 
results  of  several  computational  experiments  have  been  incorporated  into  a 
mantiscript  which  has  been  accepted  for  publication  in  ('ompul.  Me.thti.  iti  Appl. 
Mcrh.  Engr.  flOj.  Wv  have  also  extende<l  our  solution,  error  estimation,  and  local 
refinement  procedures  to  two-dimensional  parabolic  problems.  Details  of  the 
method,  which  uses  dynamic  tree  structures  to  store  data,  and  the  results  of  several 
academic  and  practical  examples  were  report<>d  in  .''liatane  .4djerid's  Ph.D.  disserta¬ 
tion  [8  . 

Th('  ai)proach  of  Adjerid  and  Flaherty  jl3j  couj>les  mesh  motion  ami  error  esti¬ 
mation  to  the  solution.  This  j)roduces  a  metlnxl  with  great  numerical  stability:  how¬ 
ever.  there  are  many  instances  where  this  e.xtra  effort  is  not  necessary.  Thus,  we  are 
considering  adaptive  local  refinement  technitiues  that  use  rectangular  and  trapezoidal 
space-time  elements  (cf.  Flaherty  and  Moore  jS.Cj  and  Dieterman  et  al.  [I4j).  These 
methods  discretize  and  -olv<'  a  jtroblem  for  one  time  ste})  using  a  finite  element- 
(ialerkin  ftrocedure.  .\t  th<'  end  of  the  time  stej).  the  discretization  error  is 
estimated,  finer  subgrids  of  spac«>-time  elements  an'  added  to  regions  of  high  error, 
and  the  probh’iii  is  recursively  solv«’d  again  on  these  rt'gions.  The  process  terminates 
and  the  int(>gration  continues  to  the  next  time  stej)  when  the  estimated  error  on  each 
grid  is  less  than  a  pri'scrilx’d  toh'rance. 

Then-  is  more  to  be  done  with  our  sj»ac<'-time  local  refinement  codes.  Flaherty 
and  graduate  stmh'iits  M.  .1.  Coyh'  and  P.  K.  Moore  are  concentrating  on  developing 
hierarchic  error  estinuition  techniques  that  an'  similar  to  those  of  .\dierid  and 
Flahe:  ..  13. IG:  for  the  metlnxl  of  lines.  The  essential  ilifficulties  an'  (i)  the  absence' 
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of  any  superconvergence  in  time  and  (ii)  the  appropriate  initial  and  boundary  condi¬ 
tions  to  apply  at  coarse-fme  mesh  interfaces.  In  addition,  Jackson  and  Flaherty  [17] 
have  extended  the  finite  element  basis  to  include  discontinuous  trial  functions.  They 
used  these  in  conjunction  with  a  moving  mesh  scheme  of  Harten  and  Hyman^  to 
develop  a  method  that  is  extremely  well  suited  to  the  solution  of  hyperbolic  systems 
with  shocks.  The  method  captures  shocks  and  contact  surfaces  as  true  discontinui¬ 
ties  without  any  artificial  diffusion  or  oscillations. 

We  have  developed  an  adaptive  finite  difference  method  for  two-dimensional 
systems  of  conservation  law's.  Preliminary  work  on  a  technique  that  used 
MacCormack's  method  on  a  moving  grid  of  quadrilateral  elements  was  reported  by 
Arney  and  Flaherty  [5,15].  This  mesh  moving  technique  used  a  clustering  algorithm 
of  Berger''  to  locate  and  group  regions  of  high  error  into  rectangles  that  isolate  spa¬ 
tially  distinct  phenomena.  An  algebraic  mesh  moving  function  was  then  used  to 
move  and  align  the  mesh  with  the  regions  of  high  error.  This  procedure  is  very  effi¬ 
cient  and.  unlike  many  other  two-dimensional  mesh  moving  techniques,  it  has  no 
problem  dependent  parameters. 

Major  D.  Arney.  a  professor  at  the  T’.S.  Military  Academy,  completed  his 
Ph.D.  studies  in  1985  under  Flaherty's  direction.  He  has  combined  the  technique  of 
|5.15]  with  a  recursive  local  refinement  strategy  and  used  it  to  create  an  adaptive  fin¬ 
ite  difference  code  for  nonlinear  two-dimensional  hyperbolic  systems  [9).  The  code 
has  been  ajtplied  to  several  difficult  compressible  flow  problems  that  contain  shocks, 
contact  stirfaces.  and  expansion^.  It  appears  to  be  able-  to  utilize  the  best  features  of 
mesh  moving  and  refinement. 

2.  Recent  Interactions 

Professors  Flaherty  and  O'Malley  and  graduate  students  lectured  and/or  visited  the 
following  conferences  and  organizations  during  1985. 

J.  E.  Flaherty  visited  Los  .41amos  National  Laboratory.  Los  .\lainos.  20-21  June. 

1985.  He  held  technical  discussions  with  Drs.  J.  M.  Hyman.  T.  Manteufel.  .1. 

Dendy.  B.  Swartz,  and  B.  Wendroff  and  lectured  on  "Moving  Finite  Difference 

and  Finite  Element  .Methods  with  Local  Refinement  for  Time  Dependent  Partial 

Differential  Etpiations," 


*  .\  Harten  and  .1  M  ll>tiian.  ’’Self  Xdiusiinp  (iri<l  Methixis  for  ( tno-Dimensional  Hyperbolic  ('on- 
servation  Laws."  J.  C\)iii|i.  Phys..  \  o|,  .Xa,  'JOP-r/J.-i.  I'js.'! 

■  M  J  Hereer.  ".Xdajitue  Mesh  Kerinenieni  for  Hyperliolir  I’arlial  Differential  Lquaiiotis."  Report 
No  .S'r.X N-( ','s-s2-U2  t.  Departnient  of  (  oiiipiiter  Sciencx,  .Stanford  I  niversity.  iyS2 
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J.  E.  Flaherty  attended  the  ASME/ASCE  Applied  Mechanics  Conference.  Albu- 
querqiie,  24-26  June.  1985.  He  presented  an  invited  lecture  on  "Local  Refine¬ 
ment  and  Moving  Finite  Element  Methods  for  Parabolic  Partial  Differential 
Equations"  as  part  of  the  .Symposium  on  Recent  .Advances  in  Computational 
Mechanics.  He  also  had  technical  discussions  with  Profs.  R.  Ewing  of  Wyoming 
University.  J.  T.  Oden  of  the  I'niversity  of  Te.xas.  and  R.  D.  Russell  of  Simon 
Frasier  University. 

R.  E.  O'Malley  attended  the  AMS  .''yrnposium  on  (.'ombustion  Theory  at  Cor¬ 
nell  University  for  a  week  in  July.  Especially  valuable  correspondence  concern¬ 
ing  singular  models  for  chemical  kinetics  resulted  with  Professor  S.-H.  Lam  of 
Princeton  University. 

R.  E.  O'Malley  visited  Los  .4iamos  National  Laboratory  and  its  Center  for  Non¬ 
linear  Studies  for  a  week  in  .August.  He  had  technical  discussions  about  singular 
perturbations  with  B.  Nicolaenko.  P.  Fife,  D.  Brown.  D.  Cohen.  J.  M.  Hyman. 
D.  .McLaughlin,  and  J.  D.  .Murray,  among  others. 

J.  E.  Flaherty  co-organized  (with  D.  A.  Drew)  the  Conference  on  Mathematics 
.\pplied  to  Fluid  Mechanics  and  Stability  which  was  dedicated  in  memory  of 
Professor  R.  C.  Di  Prinia  and  held  at  RPl.  9-11  September.  1985. 

J.  E.  Flaherty  visitied  the  IBM  T.  J.  Watson  Research  Center.  Yorktown 
Heights.  9  October.  1985  and  lecttired  on  ".\daptive  Finite  Difference  and  Finite 
ElenieiiT  Methods  for  Time  Dej)endent  Partial  Differential  Equations."  He  also 
liad  technical  discussions  with  Drs.  W.  Liniger  and  F.  Odeh. 

.1.  E.  Flaherty  visited  the  ('ourant  Institute  of  Mathematical  Sciences.  New 
York.  IS  October.  1985.  He  had  technical  discussions  with  Profs.  M.  Berger 
and  P.  Colh'la  and  lectured  on  ".XdajJtive  Finite  Difference  and  Finite  Element 
.M  «’thods  lor  Time  Dependent  Partial  Differential  Eijuations." 

S.  ,\djerid  and  .1.  E.  Flaherty  attended  the  SI,\M  Fall  meeting,  .\rizona  State 
University.  Temi)e.  28-30  Octob<T.  1985.  Thc'y  had  technical  discussions  with 
s<'verai  individuals,  including  Dr.  R.  Y.  Chin  of  Lawrence  Livermore  Labora¬ 
tories.  Prof.  R.  Bank  of  the  Univt'rsity  of  California  at  San  Diego,  and  Prof.  R. 
D.  Rusx'll  of  Simon  Frash-r  University. 

ProfesMir  O'.Malley  also  atfemh-d  the  SI.\M  Fall  met'ting  in  Tempe.  This  pro- 
vid<'d  an  o[>porf unify  for  dis^u•^sion^  about  numerical  metluxls  for  two-point  stiff 
problem'  with  (i.  Dahhpii't.  R.  Russell.  C.  Ringhofer.  and  ('.  Schmeisser. 
iimong  other'.  In  Octolier.  he  gave  a  seminar  on  the  subject  at  the  University 
of  .Montreal  anil  a  'eric'  of  thr<‘e  h’ctures  at  tin’  Uni^’ersity  of  Ne%v  Hampshin'. 


S.  Adjerid,  J.  E.  Flaherty,  and  P.  K.  Moore  attended  the  SIAM  Conference  on 
Parallel  Processing  for  Scientific  Computing,  Norfolk,  18-21  November,  1985. 
They  had  technical  discussions  with  several  individuals  including  Maj.  D.  C. 
Arney  of  the  United  States  Military  Academy.  Drs.  R.  C.  Y.  Chin  and  L.  Pet- 
zold  of  Lawrence  Livermore  Laboratories.  Dr.  S.  F.  Davis  of  the  Naval  Surface 
Weapons  Onter.  and  Prof.  J.  Van  Rosendale  of  the  University  of  Utah. 

S.  ,\djerid.  J.  E.  Flaherty,  and  P.  K.  Moore  attended  a  day-long  workship  on 
adaptive  methods  for  partial  differential  equations  at  the  University  of  Mary¬ 
land  on  24  Januaiy  1986.  Each  of  them  presented  a  lecture  on  adaptive 
methods.  .41so  present  at  the  workshop  were  I.  Babuska,  J.  Osborne,  and  E. 
Rank  of  the  University  of  Maryland:  P.  Baehman  and  M.  Shephard  of  RPI;  J. 
Thomas  of  .\FOSR:  S.  F.  Davis  and  W.  Szymezak  of  the  Naval  Surface 
Weapons  (Center;  and  other  students  and  faculty. 

Professor  O'Malley  is  chairman  of  the  1986  SL4M  National  Meeting.  Its  pro¬ 
gram  will  feature  many  sessions  relevant  to  this  contract.  He  is  also  organizing 
a  SI.WI  Conference  on  Numerical  Methods  for  Singularly  Perturbed  Problems. 

Effective  .January  1.  1986.  Professor  O’Malley  is  SIAM's  Vice  President  for  Pub¬ 
lications.  He  has  also  joined  the  editorial  board  of  the  IMA  Journal  on  .\pplied 
Mathematics. 

Drs.  O'Malley  and  Nefedov  visited  (California  in  Febniary  to  discuss  singular 
jjorturbation  problems  with  Professors  (Cohen.  Lorenz.  Lagerstrom.  and  Kreiss  at 
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ON  THE  NUMERICAL  SOLUTION  OF  SINGULARLY-PERTURBED  VECTOR 

BOUNDARY  VALUE  PROBLEMS 


.  Robert  E.  O'Malley.  Jr. 

and 

Joseph  E.  Flaherty 

Department  of  Mathematical  Sciences 
Rensselaer  Polytechnic  Institute 
Troy.  New  York  12180  US.4 

ABSTRACT 

Numerical  procedures  are  developed  for  constructing  asymptotic  solutions  of 
certain  nonlinear  singularly-perturbed  vector  two-point  boundary  value  problems 
having  boundart'  layers  at  one  or  both  end  points.  The  asymptotic  approximations 
are  generated  numerically  and  can  either  be  used  as  is  or  to  furnish  a  tw’o-point 
boundary  value  code  (e.g..  COLSYS)  with  an  initial  approximation  and  a  nonuni¬ 
form  computational  mesh.  The  procedures  are  applied  to  a  model  problem  that  indi¬ 
cates  the  possibility  of  multiple  solutions  and  problems  involving  the  deformation  of 
a  thin  nonlinear  c'lastic  beam  resting  on  a  nonlinear  elastic  foundation. 
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ADAPTIVE  FINITE  ELEMENT  METHODS 
AND  THE  NUMERICAL  SOLUTION  OF  SHEAR  BAND  PROBLEMS 


Donald  A.  Drow  and  .losrph  E.  Flaherty 
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Rensselaer  Polytechnic  Institute 
Troy.  New  York  12180 

ABSTRACT 

Shear  hands  are  localized  regions  of  very  high  shear  strain  which  arise  as  a  result  of 
high  rates  of  loading.  They  occur  in  metal  forming  and  cutting  processes  and  in 
impact  and  penetration  prohlems.  In  this  paper,  we  descril)e  a  model  for  the  forma¬ 
tion  of  shear  hands  in  simple  shear  that  involves  the  description  of  irreversible 
mechanical  shear  and  the  resulting  heat  release.  The  location  of  a  shear  hand  is  unk¬ 
nown  in  advance,  and  the  evohition  results  in  large  gradients  of  displacement,  velo¬ 
city.  and  temperature.  Shear  hand  formation,  therefore,  offers  an  interesting  and 
physically  important  application  of  a  code  able  to  resolve  small-scale  transient  struc¬ 
tures.  In  this  paper,  we  use  an  adaptive  finite  element  code  to  solve  .several  prob¬ 
lem'  involving  shear  hand  formation.  The  code  automatically  locates  regions  with 
large  gradients  and  adaptively  concentrates  finite  elements  there  in  order  to  minim¬ 
ize  ap{)roximately  the  development  of  shear  hands  uiuh'r  many  circumstances  and 
indicate  some  possible  mechanisms  for  their  formation. 
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ABSTRACT 

We  discuss  an  adaptive  finite  element  method  for  solving  initial- boundary  value 
problems  for  vector  systems  of  partial  differential  equations  in  one  spare  dimension 
and  time.  The  method  uses  piecewise  bilinear  rectangular  space-time  finite  elements. 
For  each  time  step,  the  grid  is  automatically  refined  in  regions  where  the  local 
discretization  error  is  estimated  as  being  larger  than  a  prc'scribed  tolerance.  We  dis¬ 
cuss  several  aspects  of  our  algorithm,  inchiding  the  tree-  structure  that  is  used  to 
represent  the  finite  element  solution  and  grids,  an  error  estimation  technique,  and 
initial  and  hotindaiy  conditions  at  coarse-fine  mesh  interfaces.  We  also  present  the 
results  of  several  computational  examples  an<l  ex]>eriment'. 
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ABSTRACT 

The  phenomenon  of  quenchinj'  of  cylindrical  convergen  flames  is  considered  in  the 
framework  of  constant  density  approximation.  The  parameters  dependence  of  the 
effect  is  studied  numerically  using  an  adai)tive  finite  element  method.  The  numeri¬ 
cal  results  are  in  a  good  agreement  vvitli  jiredictions  of  theoretical  analysis  from 
(Frankel  and  Sivashinsky.  1984]. 
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ABSTRACT 

Wo  discuss  an  adaptive  mesh  moviii};  terhnicjue  that  i"  used  witli  a  finite  difference 
or  finite  element  scheme  to  solve  initial-boundary  value  jirobleiii^  for  vector  system^ 
of  partial  differential  equations  in  two  space  tlimensions  and  time.  The  mesh  inoviii}; 
technique  is  based  on  an  algebraic  node  movement  function  determined  from  the 
j>ropagafion  of  significant  error  regions.  The  algoritlnn  is  (h-signed  to  b('  flexible,  so 
that  it  can  be  used  with  many  e.xisting  finite  difference  or  finite  element  methods. 
To  test  the  algorithm.  w<'  iiniilemented  the  mesh  inov<’r  in  a  system  rode  along  with 
an  initial  mesh  generator  and  a  MacCorniack  finite  volume'  integrator  to  solve  hyper¬ 
bolic  vector  systems.  Results  are  presented  for  several  coiuiiutational  examples.  The 
moving  mesh  n’duce^  disjx-rsion  errors  m-ar  shocks  and  wave-  fronts  and  thereby  can 
reduce  the  grid  re<iuirements  necessary  to  compute  accurate  solutions  while  increas¬ 
ing  computational  efficienc_\  . 
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ABSTRACT 

We  discuss  an  adaptive  finite  element  method  for  solving  initial-boundary  value 
problems  for  vector  systems  of  partial  differential  equations  in  one  space  dimension 
and  time.  The  method  uses  piecewise  bilinear  rectangular  space-time  finite  elements. 
For  each  time  step,  the  grid  is  automatically  refined  in  regions  where  the  local 
discretization  error  is  estimated  as  being  larger  than  a  prescribed  tolerance.  We  dis¬ 
cuss  several  aspects  of  our  algorithm,  including  the  tree  structure  that  is  used  to 
represent  the  finite  element  sohition  and  grids,  an  error  estimation  technique,  and 
initial  and  boundary'  conditions  at  coarse-fine  mesh  interfaces.  We  also  present  the 
results  of  several  computational  examples  and  exi»eriments. 
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ON  THE  SIMULTANEOUS  USE  OF  ASYMPTOTIC 
AND  NUMERICAL  METHODS  TO  SOLVE  NONLINEAR 
TWO-POINT  PROBLEMS  WITH  BOUNDARY  AND  INTERIOR  LAYERS 
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Troy.  New  York  12180 

ABSTRACT 

The  purpose  of  this  paper  is  to  provide  a  survey  concerning  boundary  value  problems 
for  certain  systems  of  nonlinear  singularly  perturbed  ordinary  differential  equations. 
The  aim  is  to  emphasize  important  and  difficult  oi)en  problems  needing  more  study, 
in  terms  of  both  mathematical  and  numerical  analysis  and  computational  experi¬ 
ments.  Some  currently  tractable  problems  are  discussed  in  detail. 
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ADAPTIVE  FINITE  ELEMENT  METHODS  FOR  TIME 
DEPENDENT  PARTIAL  DIFFERENTIAL  EQUATIONS 
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Department  of  Mathematical  Sciences 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  12180 

ABSTRACT 

VVe  discuss  adaptive  finite  element  methods  for  solving  initial-boundary  problems  for 
vector  systems  of  partial  differential  equations  in  one  and  two  space  dimensions. 

One  dimensional  systems  are  discretized  using  jiiecewise  linear  finite  element  approx¬ 
imations  in  space  and  a  backward  difference  code  for  stiff  ordinary  differential  sys¬ 
tems  in  time.  A  spatial  error  estimate  is  calculated  using  piecewise  quadratic 
approximations  that  use  the  superconvergence  property  of  parabolic  systems  to  gain 
computational  efficiency.  This  error  estimate  is  used  to  move  and  refine  the  finite 
element  mesh  in  order  to  equidistribute  a  measure  of  the  total  spatial  error  and  to 
satisfy  a  pre.scribed  error  tolerance.  Ordinan-  differential  equations  for  the  spatial 
error  estimate  and  the  mesh  motion  are  integrated  in  tiiiK'  using  the  same  backward 
difference  software  that  is  n.sed  to  determine  th<'  numerical  solution  of  the  partial  dif¬ 
ferential  system.  The  one-dimensional  algorithm  combines  mesh  moving  with  local 
refinement  in  a  relatively  efficient  manner  and  attemjits  to  eliminate  problem- 
dependent  numerical  parameters.  A  variety  of  examples  that  motivate  our  mesh 
moving  strategy  and  illustrate  the  performance  of  our  algorithms  are  presented. 

Two-dimensional  systems  are  discretize*!  using  jiic’cewisr  l)iiinear  finit*'  element 
approximations  in  space  and  the  .same  backward  difference  software  that  is  used  for 
one-dimensional  systems  in  time.  A  spatial  error  estimate  i><  calculated  using  piece- 
wise  cubic  approximations  that  use  the  superconvergence  jiroperty  of  parabolic  sys¬ 
tems.  This  error  estimate  is  used  to  locally  refnu'  ;»  stationary  finite  eieiiu'nt  mesh  in 
order  to  satisfy  a  prescribed  spatial  error  toleranci . 

Several  linear  and  nonlinear  examples  are  presented  to  illustrat*'  the  effectiveness  of 
our  error  estimation  technique  and  the  performanc«>  of  our  adaptive  algorithms. 
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ABSTRACT 

Wo  discuss  an  adaptive  mesh  algorithm  that  can  be  used  with  a  finite  difference  or 
finite  element  scheme  to  solve  initial-boundary  value  problems  for  vector  systems  of 
time  dependent  partial  differential  equations  in  two  space  dimensions.  Our  algo¬ 
rithm  combines  the  adaptive  techniques  of  mesh  moving,  static  rezoning,  and  local 
mesh  refinement.  The  nodes  of  a  coarse  mesh  of  quadrilateral  cells  are  moved  by  a 
simple  algebraic  node  movement  function,  determined  from  the  geometry  and  propa¬ 
gation  of  regions  having  statically  significant  discretization  error  or  mesh  movement 
indicators.  The  local  mesh  refinement  method  recursively  divides  cells  of  the  moving 
coarse  mesh  within  clustered  regions  that  contain  nodes  with  large  error  until  a  user 
prescriln'd  error  tolerance  is  satisfied.  These  finer  grids  are  properly  nested  within 
the  moving  coarse  mesh  to  provide  for  simpler  «lata  structures  and  interface  condi¬ 
tions  between  the  fine  and  coarse  meshes. 

Our  {iroceclun’  is  designed  to  h(>  tlexible,  so  that  it  can  be  used  with  many  existing 
finite  difference  and  finite  element  schemes  and  with  different  error  estimation  pro¬ 
cedures.  To  test  our  adaptive  mesh  algorithm,  we  implemented  it  in  a  system  code 
with  an  initial  mesh  generator,  a  MacCormack  finite  difference  scheme  for  hyperbolic 
vector  syst<>nis  of  conservation  laws,  ami  a  Richard.son  extrapolation  based  error  esti¬ 
mation.  Results  are  presented  for  several  computational  e.xaniples. 

The  moving  mesh  technique  reduces  disfiersive  errors  near  shocks  and  wave  fronts. 
Therefore,  it  reduces  the  grid  retpiirements  necessary  to  compute  accurate  solutions 
and  thus  increases  computational  efficiency.  The  local  mesh  refinement  provides 
'iiialler  mesh  spacings  and  time  steps  in  regions  where  the  problem  is  difi'iciilt  to 
solve,  thus  providing  increased  accuracy  and  enabling  t'rror  tolerances  to  be 
achieved. 

I’h  I)  Dis^rrtatiofi.  Dcparlmcril  nf  .Vlatlicrrialical  .sciences.  Hens.selaer  I’olyechnic  In.slitiiie.  Tr<)\ 


DECOUPLING  AND  BOUNDARY  VALUE 
PROBLEMS  FOR  TWO-TIME-SCALE  SYSTEMS 

R.  E.  O'Malley.  Jr. 

Department  of  Mathematical  Sciences 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  12180 
and 

R.  M.  M.  Mattheij 
Catholic  University 
Nijmegen,  The  Netherlands 

ABSTRACT 

This  paper  presents  an  analytical  technique  for  decoupling  slow  and  fast  dynamics 
for  linear  systems  of  ordinary  differential  equations.  The  method,  based  on  Riccati 
transformations,  is  of  substantial  computational  value  for  solving  stiff  two-point 
boundary  value  problems. 
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ABSTRACT 

We  study  the  stability  of  several  mesh  equidistribution  schemes  for  time  dependent 
partial  differential  equations  in  one  space  dimension.  The  schemes  move  a  finite 
difl'erence  or  finite  element  mesh  so  that  a  given  quantity  is  uniform  over  the 
domain.  We  consider  mesh  moving  methods  that  are  based  on  solving  a  system  of 
ordinary  differential  equations  for  the  mesh  velocities  and  show  that  many  of  these 
methods  are  unstable  with  respect  to  an  equidistributing  mesh  when  the  partial  dif¬ 
ferential  system  is  dissipative.  Using  linear  perturbation  techniques,  we  are  able  to 
develop  simple  criteria  for  determining  the  stability  of  a  particular  method  and  show 
how  to  construct  stable  differential  systems  for  the  mesh  velocities.  Several  exam¬ 
ples  illustrating  stable  and  unstable  mesh  motions  are  presented. 
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ABSTRACT 

In  two  recent  papers.  Slemrod  has  suggested  that  the  well  known  Lax-Friedrichs 
finite  difference  method  may  provide  a  natural  method  for  the  numerical  integration 
of  initial  value  problems  with  an  anomalous  equation  of  state,  e.g..  a  van  der  Waals 
fluid.  In  this  note  we  review  these  ideas  and  present  the  results  of  a  numerical 
experiment  which  attempts  to  simulate  the  dynamics  of  a  van  der  Waals  like  fluid. 
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Dedicated  in  Memory  of  Richard  C.  DiPrima 
ABSTRACT 

We  discuss  a  moving  finite  element  method  for  solving  vector  systems  of  time  depen¬ 
dent  partial  differential  equations  in  one  space  dimension.  The  mesh  is  moved  so  as 
to  equidistribute  the  spatial  component  of  the  discretization  error  in  Wo 

present  a  method  of  estimating  this  error  by  using  p-hierarchic  finite  elements.  The 
error  estimate  is  also  used  in  an  adaptive  mesh  refinement  procedure  to  give  an  algo¬ 
rithm  that  combines  mesh  movement  and  refinement.  We  di.screTize  the  partial  dif¬ 
ferential  equations  in  space  using  a  Galerkin  procedure  with  piecewise  linear  ele¬ 
ments  to  approximate  the  solution  and  quadratic  elements  to  estimate  the  error.  .4 
system  of  ordinan.'  differential  equations  for  mesh  v«'lociti('s  are  used  to  control  ele¬ 
ment  motions.  We  use  existing  software  for  stiff  ordinarj'  differential  equations  for 
the  tpmi)oral  integration  of  the  solution,  the  error  estimate',  and  the  mesh  motion. 
(h)iiiputational  results  using  a  code  based  on  our  methoel  are  presented  for  several 
examples. 
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ABSTRACT 

We  consider  two  adaptive  finite  element  techniques  for  iiarabolic  partial  differential 
equations  (PDEs)  that  are  based  on  using  error  estimates  to  control  mesh  refine¬ 
ment.  One  technique  is  a  method  of  lines  (MOLl  apjiroach  that  uses  a  (talerkin 
method  to  di.scretize  the  PDEs  in  sjiace  and  inijilicit  multi-stej)  integration  in  time. 
Spatial  elements  are  added  and  deleted  in  regions  of  high  and  low  error  and  are  all 
advanced  with  the  same  sequence  of  varying  time  steps.  The  second  technique  is  a 
local  refinement  method  (LRM)  that  uses  (talerkin  approximations  in  both  space 
and  time.  Fine  grids  of  space-time  elements  are  added  to  coarser  grids  and  the  prob¬ 
lem  is  recursively  solved  in  regions  of  high  error. 
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ABSTRACT 

We  discuss  an  adaptive  mesh  moving  technique  that  can  be  tised  with  a  finite  ele¬ 
ment  difference  or  finite  element  scheme  to  solve  initial- boundary  value  problems  for 
vector  systems  of  partial  differential  equations  in  two  <pace  dimensions  and  time. 
The  mesh  moving  technique  is  based  on  an  algebraic  node  movement  function  deter¬ 
mined  from  the  geometry  and  propagation  of  regions  having  significant  discretization 
error  indicators.  Our  procedure  is  designed  to  be  flexible,  so  that  it  can  be  used  with 
many  existing  finite  difference  and  finite  element  methods.  To  test  the  mesh  moving 
algorithm,  we  implemented  it  in  a  system  code  with  an  iidtial  mesh  generator  and  a 
.MacC'ormack  finite  difference  scheme  on  quadrilateral  cells  tor  hyperbolic  vector  sys¬ 
tems  of  conservation  laws.  Results  are  presented  for  si-veral  computational  exam¬ 
ples.  The  moving  mesh  scheme  reduces  dispersive  errors  near  shocks  and  wave 
fronts  and  thereby  reduces  the  grid  requirements  nect'ssary  to  compute  accurate  solu¬ 
tions  while  increasing  computational  efficiency. 
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ABSTRACT 

We  discuss  a  moving  mesh  finite  element  method  for  solving  initial-boundary  value 
problems  for  vector  systems  of  partial  differential  equations  in  one  space  dimension 
and  time.  The  system  is  discretized  using  piecewise  linear  finite  element  approxima¬ 
tions  in  space  and  a  backward  difference  code  for  stiff  ordinary  differential  systems  in 
rime.  spatial  error  estimation  is  calculated  using  piecewise  quadratic  approxima¬ 
tions  that  use  the  superconvergcnce  properties  of  parabolic  systems  to  gain  computa¬ 
tional  efficiency.  The  spatial  error  estimate  is  used  to  move  and  locally  refine  the 
finite  element  mesh  in  order  to  equidistribute  a  measure  of  the  total  spatial  error  and 
To  satisfy  a  prescribed  error  tolerance.  Ordinary  differential  equations  for  the  spatial 
error  estimate  and  the  me>ii  motion  are  integrated  in  time  using  the  same  backward 
diff('renc(>  software  that  is  used  to  determine  the  numerical  solution  of  the  partial  dif¬ 
ferential  system. 

We  present  several  details  of  an  algorithm  that  may  be  used  to  develop  a  general 
purpose  finite  element  code  for  one-dimensional  parabolic  partial  differential  systems. 
The  algorithm  combines  niesli  motion  aiul  local  refineuKuit  in  a  relatively  efficient 
manner  and  attempts  to  idiminate  problem-dejxuident  nuriK'rical  parameters.  A 
variety  of  examples  that  motivate  our  mesh  moving  strategy  and  illustrate  the  per¬ 
formance  of  our  algorithm  are  presented. 
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abstract 

We  develop  an  adaptive  finite  element  method  for  solving  nonlinear  hyperbolic  sys¬ 
tems  of  conservation  laws  in  one  space  dimension  and  time.  The  method  uses 
discontinuous  piecewise  linear  trial  functions  and  continuous  piecewise  cubic  test 
functions  on  a  moving  mesh  of  triangular  space-time  elements.  The  mesh  is  moved 
by  a  technique  that  uses  a  weighted  average  of  the  local  characteristic  speeds  to 
select  nodal  velocities.  We  show  that  discontinuous  finite  element  approximations 
on  a  moving  mesh  have  the  potential  of  accurately  resolving  physical  discontinuities 
in  the  solution.  Several  computational  examples  are  presented  to  illustrate  the  per¬ 
formance  of  the  method. 
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ABSTRACT 

The  study  of  singularly  perturbed  initial  value  problems  for  nonlinear  systems  of 
ordinary  differential  equations  parallels  the  analysis  underlying  the  development  of 
numerical  algorithms  for  obtaining  solutions  of  such  systems  of  stiff  differential  equa¬ 
tions.  This  paper  seeks  to  emphasize  the  advantages  of  combining  these  two  sub¬ 
stantial  research  efforts.  It  develops  insight  and  intuition  based  on  a  sequence  of 
solvable  model  problems,  and  it  relates  a  variety  of  literature  scattered  throughout 
publications  concerned  with  asymptotic  and  numerical  analyses,  stability  and  control 
theory,  and  specific  topics  in  applied  mathematical  modeling. 
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